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SUMMARY 


An experimental and analytical study has been made of some features 
of the turbulent heat diffusion behind a line heated wire stretched 
perpendicular to a flowing isotropic turbulence. The mean temperature 
distributions have been measured with systematic variations in wind 
speed, size of turbulence-producing grid, and downstream location of 
heat source. The nature of the temperature fluctuation field has been 
studied. 

A comparison of Lagrangian and Eulerian analyses for diffusion in 
a nondecaying turbulence yields an expression for turbulent-heat-transfer 
coefficient in terms of turbulence velocity and a Lagrangian "scale.'.' 

The ratio of Eulerian to Lagrangian microscale has been determined 
theoretically by generalization of a result of Heisenberg and, with 
arbitrary constants taken from independent sources, shows rough agree- 
ment with experimental results. 

A convenient form has been deduced for the criterion of inter- 
changeability of instantaneous space and time derivatives in a flowing 
turbulence. 


INTRODUCTION 


One of the most striking aspects of turbulent motion in fluids is 
its dispersive property. This "convective diffusion," illustrated by 
the general statistical tendency of (noncontiguous) fluid elements to 
get farther apart with increasing time, was probably first observed long 
before the era of analytical fluid mechanics. An analytical start on 
this problem was not made, however, until the now-classic work by 
Taylor in 1921 on diffusion by continuous movements (reference 1) . Not 
only did this paper lay a groundwork for the study of turbulent diffusion 
but it also represented a forward step in the ideas essential to develop- 
ment of a general statistical theory of turbulence, a field which had . 
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scarcely progressed since Reynolds' original- formulation of the equations 
of motion for a flow in which mean and fluctuating parts could be 
distinguished. 

The diffusive action of a turbulent flow may manifest itself in 
various ways, depending upon the initial and/or boundary conditions and 
upon the interests of the observer. The following possible measures of 
the diffusive powers are neither exhaustive nor mutually independent: 

(1) The average rate of dispersion of particles from a fixed source 

(2) The average rate of increase of spacing between different 

particles 

(3) The average rate of transport of particle concentration under 

a given mean concentration gradient 

(4) The average rate of increase of the length of a fluid line 

(5) The average rate of increase of the area of a fluid surface 

The word "particles" means simply indelibly tagged fluid elements, much 
smaller than the smallest length associated with the turbulence. 

The present report is concerned primarily with measure (l). The 
measurements have all been made in the thermal wake of a long thin heated 
wire mounted perpendicular to an isotropic turbulent air flow and 
producing no turbulent wake. Here the tagging is thermal, and the degree 
of indelibility (negligibility of molecular diffusion) is one of the 
matters to be investigated. 

The diffusive property (for scalars) of a turbulent flow is appar- 
ently a secondary characteristic at least in the sense that it need not 
explicitly enter the dynamical problem. The diffusion may be regarded 
as a kinematic phenomenon, to be deduced from the dynamical solution to 
the problem if and when the latter is obtained. Thus the objective of 
research on turbulent diffusion may be to seek a connection between the 
diffusive and the dynamical statistical variables, even before the 
complete dynamical theory is available. 

Measure ( 3 ) is usually termed the "turbulent transport" or "transfer" 
problem. Although of extensive practical importance, it has not yet been 
subjected to genuine theoretical study. 

Most of the semiempirical "theories" of turbulent transport, for 
both scalar and vector properties, employ an Eulerian formulation of the 
basic equations, and up to now they have been unable to relate the 
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turbulent transport correlation to other statistical functions describing 
the flow. Taylor (reference l) showed that in the simple case of a 
homogeneous field of isotropic turbulence, and even in a decaying 
isotropic turbulence (reference 2), a Lagrangian formulation of the 
transport (i.e., diffusion} the terms will be used interchangeably) 
problem leads to some important results. ^ 

Up to the present time little theoretical or experimental work has 
been done to find relations, if any, between the Lagrangian statistical 
measures of a turbulent field with its Eulerian statistical measures. 

Since turbulence dynamics seems best handled in the latter terms and 
turbulent diffusion in the former, it is evident that such a connection 
is important. Hence, one of the purposes of the present experiments has 
been to compare the magnitudes of some of these quantities under varia- 
tions in the turbulent field. For example, the postulates of Taylor 
and Heisenberg on a relation between Lagrangian and Eulerian microscales 
can be examined and, in corrected form, compared with experiment. 

The turbulent diffusion from a fixed line source can be set up 
analytically as an ordinary (Eulerian) "heat-transfer" problem, permitting 
a start to be made in relating measures (l) and (3) of the diffusive 
power of a turbulent flow, under certain simplifying assumptions. 

Measures (4) and (especially) ( 5 ) may well be classed as character- 
istic of the "turbulent -mixing" problem rather than diffusion in the 
common connotation. 

Experimental work on diffusion from a fixed local source in a 
turbulent flow has been meager. In isotropic turbulence, there have 
been the measurements of Schubauer (reference 4 ), Simmons (reported by 
Taylor in reference 2), Dupuis (reported by Kampe de Feriet in refer- 
ence 5), Frenkiel (reference 6 ), and Collis (reference 7 )* Of these, 
only the data of Simmons and Collis are extensive enough to permit 
confident computation of the Lagrangian correlation function. In 
turbulent shear flow, Skramstad and Schubauer (reference 8), Dryden 
(reference 9 ), and the present authors (reference 10) have measured 
distributions close to a source; Kalinske and Pien (reference 11) and 
Van Driest (reference 12) have made measurements somewhat farther 
downstream. 

None of these studies was repeated with a systematic variation of 
the properties of the turbulence. In spite of the poor precision 


lln his tensorial generalization of Taylor f s work on case (l), 
Batchelor (reference 3 ) Has chosen to call this an "Eulerian" analysis, 
describing only case (2) as "Lagrangian." In keeping with previously 
accepted nomenclature, both cases are Lagrangian, (l) involving a single 
particle, and (2) involving a pair. In fact, case (2) might be termed 
a mixed (Eulerian and Lagrangian) problem. 
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inherent in this type of measurement, it was hoped that such an approach 
would at least show up some general trends in the relations between 
Eulerian and Lagrangian variables. 

This investigation has been conducted at the Aeronautics Department 
of the John Hopkins University under the sponsorship and with the 
financial assistance of the National Advisory Committee for Aeronautics. 
The authors would like to acknowledge the assistance of Messrs. Alan 
Kistler, George Stierhoff, and Allen Gates and Miss Patricia O'Brien, 
as well as the helpful criticism of Dr. Francis H. Clauser and Dr. C. C. 
Lin. 


SYMBOLS 


M grid mesh size 

x distance downstream from grid 

x 0 . location of heating wire 

Ax = x - x Q 

y 

z 
r 
U 
V 
u 
v 
w 

u* = {W- 

V T = 

c root -mean- square molecular velocity 

instantaneous temperature (measured above ambient room 
temperature) •. 


distance in direction of measured diffusion 
distance in direction of heating wire 
scalar distance between two points 
mean velocity in x-direction 
mean velocity in y-direction 

instantaneous velocity fluctuation in x-direction 
instantaneous velocity fluctuation in y-direction 
instantaneous velocity fluctuation in z -direction 


0 
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© 

©o 


■d 

*0 


mean temperature 

maximum mean temperature at a cross section, a function 
of Ax 

instantaneous temperature difference 
temperature difference of rectangular heat pulse 


0’ = 

j 

S 


width of rectangular heat pulse 
average on-center spacing of pulses 


n(t 0 ,t) 



v'( I) d£ 


t 


time 


T 


time difference 


f(r) Eulerian velocity correlation coefficient (notation of 

Von Karraan and Howarth) 

Ry( T ) Lagrangian correlation coefficient for nondecaying 

turbulenc e 


t^t-T = 


v(t)v(t - t) 
v’(t)v’(t - t) 


Rt]( t i) 


Lagrangian correlation as a function of for 
nondecaying and decaying turbulence 


R 



Y 


X 


Y* 



distance traveled in y-direction by a fluid particle 

distance traveled in x-direction by a fluid particle 

root -mean- square displacement of a fluid particle in 
y-direction 
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X 

Xy 


X 


*1 


Ll 

L 


A 

l 

C P 

P 

V 

k 

P 

K 


I = 


nx) 


root -mean- square displacement of a molecule 


Eulerian microscale of turbulence (k = 


1 


rf"(0) y 


Lagrangian microscale of turbulence for nondecaying 


turbulence \X y = 


1/21 


-V(o) 


Lagrangian microscale of turbulence for nondecay Jig and 


and decaying turbulence = 






Lagrangian scale IL^ = R v (t ) dT j 


JO 


Eulerian scale (L = f(r) dr 

Jo ) 


Lagrangian scale for nondecaying and decaying turbulence 

dnj 

mean free path of a molecule 
mixing length 

specific heat per unit volume at constant pressure 
density 

kinematic viscosity 
thermal conductivity 
static pressure 

dimensionless empirical constant 
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poo 


H = 2pc p U 


© dy 


H*s 


kr|i 


H 


2pc T 


00 

f(5) 

i A 


d| 


JO 


turbulent-heat-transfer coefficient 


Mi = 


T) R-qCr)) dtj 



p(d) probability density of fluctuating temperature 

Rl turbulence Reynolds number based on Eulerian scale 



Ri 


n 

max 

min 


turbulence Reynolds number based on Eulerian microscale 



mean value or ensemble average 
maximum, used as subscript 
minimum, used as subscript 
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EQUIPMENT AND PROCEDURE 
Aerodynamic Equipment 


The wind tunnel (fig. l) is an open-return NPL type tunnel with 
a 2- by 2-foot working section and a free-stream turbulence level of 
v'/U = 0.06 percent and u'/lJ = 0.05 percent at a mean velocity of 
26 feet per second. The turbulence-producing grids were as follows: 


Designation 

Type 

Mesh size 
(in.) 

Rod 

diameter 

(in.) 

Solidity 

1-in. grid 
l/2-in. grid 
l/4-in. grid 

Biplane , wood 
Woven, metal ’ 
Woven, metal 

1.00 

a .50 by O .53 
.25 

0.25 

.120 

.060 

0.437 
.41 
.42 . 


^he 0.50-in. mesh was set in direction of measured diffusion. 


They were mounted in turn at the upstream end of the working section. 

The heat source was an 0.008-inch-diameter platinum wire stretched 
vertically across the tunnel at various distances from the grid. It was 
heated by direct current to temperatures between 500° and 700° C, with 
the latter figure only at the highest operating velocity of 38.0 feet 
per second. The wire Reynolds numbers at this condition and at the two 
other velocities were as follows: 


Mean speed 
(fps) 

Wire temperature 
(°C) 

Reynolds number based on - 

Air temperature 

Wire temperature 

8.5 

500 

35 

6.4 

25.6 

500 

105 

19 

38.0 

700 

157 

19 


A preliminary investigation was made without grids to insure that 
these operating conditions did not generate a vortex street downstream 
of the heated wire. 

With the grids in place, the mean momentum wake became practically 
undetectable with total -head tube and manometer at distances greater 
than 1 or 2 inches downstream of the heating wire. 
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Measuring Equipment and Procedures 

The mean-temperature distributions were measured with a Chromel- 
Alumel thermocouple and a Leeds and Northrup type K-2 potentiometer. 

The cold junction was kept outside of the wind tunnel. 

The shadowgraph technique was used to photograph the -laminar thermal 
wake close to the source with no grid in the wind tunnel. This informa- 
tion was. applied to the problem of "correcting" the thermal wake in 
turbulent flow for the effects of molecular diffusion and of finite 
source size. Figure 2(a) is a shadowgraph of the wire wake with no grid 
in the tunnel; figure 2(b) is a typical time exposure with grid-produced 
turbulence. A resistance-thermometer traverse of the laminar wake in a 
flow of very small turbulence showed that the temperature profile had 
already become very nearly Gaussian at a distance of 1 inch (125 wire 
diameters) downstream. The white lines on the sides of the dark wake 
shadow in figure 2(a) correspond to the minimums in the second derivative 
of the density profile. For small temperature differences these coincide 
with the maximums in the second derivative of the temperature profile. 
Although the temperature differences are not small in the imm ediate 
vicinity of the wire, this condition is reasonably well satisfied at 
relatively small values of Ax as evidenced by the parabolic spread of 
this laminar wake. 

The standard deviation of the wake could then be computed from the 
spacing of these two bright lines, with the assumption of a Gaussian 
distribution. Since the closest points of traverse in the turbulent 
cases were l/2 inch (63 wire diameters) from the heat source, this was 
probably a reasonable assumption. 

As pointed out by Taylor (reference 2), the molecular and the 
turbulent diffusive phenomena are statistically independent, so that the 
squares of the standard deviations due to these two effects are additive. 
Hence the wake spread due to turbulence alone, from a true line source, 
was obtained by subtracting the square of the standard deviation of the 
laminar wake (computed from the shadowgraph) from the square of the 
standard deviation of the. total wake (from thermocouple traverses with 
grid turbulence present) at all stations. This difference was the 
square of the standard, deviation of the desired phenomena. All wake- 
spread data presented in the next section have been corrected in this 
fashion. 

Parenthetically, it should be remarked that the laminar wake in 
figure 2(a) spreads parabolically within the limits of precision, from 
at least 1 inch on, so that the effects of density differences on the 
flow phenomena must have been negligible for this investigation. 
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The transverse turbulence levels v ! /u behind the grids were 
obtained from the initial rate of spread of the mean thermal wake 
(method due to Schubauer (reference 4)) after the effects of molecular 
spread and finite source had been removed. The resulting levels were 
somewhat higher than those obtained with a hot-wire anemometer but were 
used because of their consistency with the rest of the measured dif- 
fusion curve. 

Free-stream velocity fluctuations (without grids) and the wake 
temperature fluctuations (with grids) were measured with the hot-wire 
anemometry equipment described in reference 10. The wires were 
0.00025-inch platinum etched from Wollaston; the compensated response 
of the system was flat within *2 percent over a frequency range from 
3 to 12,000 cycles per second. 

Oscillograms of the temperature fluctuations were recorded by 
photographing a blue oscilloscope tube with fast 35-millimeter film 
in a General Radio type 651-AE camera. 

Probability densities of the temperature fluctuations at fixed 
points in the mean thermal wake were determined from photodensitometer 
traverses of time-exposure photographs of a short-persistence (0.001 sec 
blue oscilloscope tube with the temperature fluctuations on one pair of 
plates and a 30,000-cycle-per-second sweep on the opposite plates. The 
technique is essentially that used by Simmons and Salter (reference 13). 


EXPERIMENTAL RESULTS 
Mean Thermal Wake . 


Complete mean-temperature wakes behind a line source of heat were 
measured for 10 different conditions. Arranged to indicate the 
systematic variation of one parameter at a time, these conditions were 
as follows: 


l-in. grid; ^ = 43.4; U, fps, 

8.5 

25-6 

38.0 

U * 25.6 fps> « 43.4; grid, in., 

1 

1/2 

1/4 

U = 25.6 fpe; ^ = 86.1} grid, in., • • 

1 

1/2 

1/4 

U = 25-6 fps} l/2-in. grid} Xq/m . . . ' 

43.4 

86.1 

172.3 

U = 25-6 fps; l/4-in. grid} xq/M 

43-4 

86.1 

172.3 
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Here U is the mean velocity, M is the grid mesh size, and x Q is 
the heat-source location measured from the grid. Since some individual 
cases enter as elements in two sequences, the total number of elements 
is more than 10. 

Two of the many mean-temperature traverses in the y-direction 
(perpendicular to mean flow and to source line) are shown in figure 3 
to give an idea of the amount of experimental scatter. The upper 
traverse was the worst of the lot, even showing an apparent skewness 
which was not borne out by the investigation as a whole. The lower 
traverse is more nearly typical of the measured temperature distribu- 
tions from which the standard deviations of the mean thermal wake were 
computed. By comparison with the reference curve, it is seen to be 
essentially Gaussian. This was the case for temperature profiles at 
all stations. Since the virtually Gaussian character of such a wake 
has already been established by several of the earlier publications, 
there seemed to be no point in reproducing here all of the large number 
of traverses measured. 

The mean-thermal-wake spread for the 10 different configurations 
studied is given in figures ^ to 7 as plots of corrected standard 
deviation Y 1 against distance from the heat source Ax. Each point 
in these figures corresponds to a complete transverse temperature 
/traverse. 

In order to have values of transverse turbulence level v'/u 
consistent with the thermal -wake behavior, these values were- determined 
from the initial angle of spread of the corrected wake standard devia- 
tion (references 2 and 4) instead of from direct hot-wire anemometer 
measurements. The results are plotted in figure 8. Since these 
represent an insufficient number of points per grid to permit the 
drawing of reliable curves, some simplifying assumptions were made based 
upon the results of several more-detailed turbulence-decay investigations. 
(See references ik to 20. The pertinent results of references lk to 1 6 
are summarized in reference 17 *) In the light of these papers, it was 
assumed that the decay curves had a common apparent origin, and that 
this was obtainable by drawing the best straight line for all the 
available points, independent of the differing wind speeds and mesh 
sizes. In computing any individual Lagrangian correlation function 
from the corresponding wake history, the turbulence decay rate was 
assumed to be given by the line drawn through this common origin and 
the specific turbulence value giving the measured initial spread angle 
for this wake. This is, of course, a very rough procedure, but the 
experimental scatter in this whole method of determining Lagrangian 
correlation functions is so great that a more extensive study of decay 
(including the resolution of inconsistencies between wake method and 
hot-wire method) seemed unwarranted at this time. 
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Temperature Fluctuations 

Distributions of temperature-fluctuation level tf'/© in the thermal 
wake have been measured by using the hot-vire anemometer as a resistance 
thermometer, that is, at a current low enough to render the sensitivity 
to velocity fluctuations negligible compared with the sensitivity to 
temperature fluctuations (reference 21). A representative distribution 
of -a'/e in the x-direction is given in figure 9* Typical transverse 
distributions are given in figures 10(a) and 10(b) . It is clear that 
the temperature-fluctuation intensity changes very little with increasing 
values of Ax. A rough explanation of the very high values of ■&'/© 
(compared with the concomitant turbulence level, for example) in terms 
of the highly intermittent structure of the thermal wake has been- given 
in reference 10 and will be discussed in more detail later in this 
report. This intermittency is shown very clearly in figure 11, a 
series of temperature oscillograms recorded at two different positions 
across the thermal wake for a fixed value of Ax and at two different 
values of Ax with y = 0. The one-sided and pulse character of the 
instantaneous temperature at a fixed point in space is also demonstrated 
by its probability density. 


THEORETICAL CONSIDERATIONS 
Homogeneous Steady Turbulence at Rest 


For a nondecaying incompressible turbulence with no mean motion, 
Taylor (reference 1) was followed in getting an expression for the mean 
time rate of diffusion in the y-direction (say) from a fixed source as 
measured by the second moment of the probability density of the diffu- 
sion, that is, the mean-square particle displacement Y^(t): 




= 2v(t) 



The bar denotes ensemble average. 

Taking v(t) inside the integral, interchanging the processes of 
integration and averaging, and introducing the Lagrangian (auto) 


) 
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correlation coefficient (following the fluid particle), where 
T = t - tp 


R V ( T ) 


v(t)v(t - t) 


there results 


Y 2 = 2V 2 



dT dT 


( 1 ) 


This is Taylor's form. Integration by parts yields a form like that in 
the work of Kampe' de Feriet (reference 5): 


( 2 ) 


In this Lagrangian analysis, v(t) is the velocity of a fluid 
particle in the y-direction at time tj v(t - t) is the velocity of 
the same particle at time t - t. Corresponding expressions can be 
written for the rate of diffusion in any direction. 

Diffusion from an infinite line source, the case to be discussed 
here, is a two-dimensional problem in the mean, and, in addition to 
equation (2), 


(t - t)R u (t) dr (3) 


A tensorial generalization of these concepts has been given by 
Batchelor (reference 3) dealing with the behavior of XjXj(t) where 
and Xj are any two of the orthogonal displacements of the particle at 
time t . 




Y 2 = 2v 2 (t - t)R v (t) dT 


It should be noted that this analysis gives no information on the 
shape of the probability density of Y(t) or of X(t). In fact there 
still exists no theory for these. However, experiments in flowing 
turbulence (the case to be considered next) show Gaussian density, 
within the experimental precision, for Y(t) at all values of t. 
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For this stationary random process, Taylor (reference 1) introduced 
the concept of the Lagrangian "scale," 


~-oo 

L l = R V (T) dr 

Jo 


( 4 ) 


These have the dimensions of time and are characteristic constants of 
the system. 

In his later work on the (Eulerian) dynamics (reference 2), Taylor 
had occasion to introduce another measure of the correlation function, 
which he called the "microscale." Applying the same geometrical concept 
to the present function, the Lagrangian microscale 


= 


•2 

rTW 


( 5 ) 


is simply the r -intercept of the vertex-osculating parabola of the even 
function R v (t). The kinematic significance is clearly shown by a 
series expansion of v(t + t) in R v (t): 


R v (t) = v(t)v(t , + ,_L ) 
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From equations (6) and (5), the Lagrangian microscale for v(t) is 


or 



(7) 


(7a) 


Turbulence in a Flowing Medium 

The dictates of both practical interest and experimental feasibility 
require analysis of the diffusion when there is a mean velocity U 
relative to the source. Since the diffusion phenomenon is linear, the 
probability density (mean- concentration distribution of tagged particles 
in the wake) is simply proportional to the superimposed probability 
densities of a continuous line of sources moving with the mean velocity 
U with their time (and space) origin at the actual fixed source. This 

is illustrated in figure 12(a) for - *. u < 1. The circles (corresponding 

to isotropy) are the standard deviations of the dispersions that would 
occur from moving sources. The envelope of these circles gives a 
measure of the mean wake. It is obvious that in general the functional 
form of the mean-concentration distribution along a line Ax = Constant 
will not be the same as the functional form of the same quantity for 
the individual source at time t, that is, at position x = Ut. 

d Y 1 

However, it is extremely likely that if « 1 (fig* 12(b)) 

the mean concentration along a line Ax = Constant becomes very nearly 
of the teame functional form as that for the source at x = Ut. For a 
zero-correlation Gaussian density, the equivalence is easily 
demonstrable . ^ 


2 The most general mathematical restrictions under which the super- 
position of a line of identical densities will yield a "cross-section” 
density of the same form has not been studied here. It is obvious that 
the condition of statistical independence is sufficient. The fact that 
the density of each of the velocity components in isotropic turbulence 
has been found to be Gaussian within the experimental precision seems 
to show that the equivalence under discussion is at least a good 
approximation. 
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The condition « 1 will always occur at large enough values 

V 1 

of Ax when — < 1- This follows from the asymptotic parabolic 
behavior, Y'(t) ~ Vt (reference 1), of the diffusive process. It will 

V 1 

occur for all* values of Ax when — « 1. Therefore, in this case a 

U 

simple approximate space-time transformation in the mean is permissible, 
and the t-variation in Taylor's theory of diffusion by continuous move- 
ments becomes a variation of x/TJ. (it must be emphasized that the 
foregoing discussion does not apply directly to the possibility of 
applying a space-time transformation to the instantaneous turbulence 
variables. This latter question will be discussed later.) 

V 1 

In the present measurements — « 1, and, therefore, the 

Ax -variation of diffusion gives an approximate measure of the Lagrangian 
correlation coefficient (in time). Equation (2) can be rewritten as; 


<e) 


Decaying Isotropic Turbulence 

When the turbulence is decaying in time (similar to space in the 

flowing turbulence) v 2 is no longer constant, and the analysis cannot 
be carried out as far as equation (l) . The same approach stops with 

^ Y 2 = v'(t) J' v* ( t - t) t R t _ T dT (9) 

where the prime denotes root -mean- square value, and 


t R t-T 


v(t)v(t - t) 
v' (t)v'(t - t) 


There is no a priori reason to believe that t R t-T is a function of t 
alone, as in the nondecaying turbulent flow. 
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At this point Taylor (reference 2) invokes the empirical fact that 
over a wide range of mean_velocities (all of which give essentially the 
same distribution of v'/U in x behind a grid) the thermal wake 
behind a line heat source at fixed x .appears to be unchanged in form, 
within the experimental error. This is consistent with dependence of 
the diffusive process upon a variable of type 


"(M 


( x i) 


Therefore, Taylor has postulated the unique dependence of -t^t-T on 

the variable T). With this postulate and the space-time transformation 
valid for small turbulence level, he arrives at 


5 l2 = 2 U 


f Y 2 = 2 

dTj 


R ti( t >i) dT ll 


where 


drj = — dx 
U 


= v* dt 


Equations (ll) and (12) look like the equation for nondecaying 
turbulence. They also give 


p 

Y = 2 


(n ~ ^R^l) ^ 1 
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A physical significance of the length r| is underscored by the 
. limiting form of equation (13) as t| — *>0 and »1. Then, 


Y' = T) 

so that T) is a measure of the lateral diffusion that would occur if 
the lateral velocity fluctuation following a particle v(t) remained 
perfectly -correlated but decreased in magnitude according to the decay 
rate of the turbulence level. 

In postulating tPt-T ~ Taylor was apparently comparing only 

diffusive processes in turbulence fields with identical ^(x). Of more 

U 

general interest is .the comparison of diffusion in fields with differing 
turbulence-level distributions. Although such a generalized application 
of his postulate is doubtless not too well applicable, it is conceivable 
it might have approximate success in the more general comparison. For 
both convenience and lack of any obviously superior alternative, his 
suggestion is therefore applied in computing the results of the measure- 
ments reported here. 3 

With the T| -postulate, the Lagrangian correlation function can 

be obtained from measurements of as a function of Ax: 


y^i) 


1 y 2 

2 dy 2 


= 1 uf. _d£_ y 2 

2 dx 2 


(14) 


3 After this work was completed. Dr. Batchelor suggested an alternative 
approximate approach: In order to construct a stationary random function 

out of the nonstationary v(t) one first normalizes the dependent 
variable with its root-mean-square value. (This has been automatically 
accomplished by use of the correlation coefficient.) If the decaying 
quantity is assumed to maintain complete similarity during decay, all 
characteristic times (e.g., Lagrangian time scale and microscale) vary 
in the same way with t and the new independent variable is constructed 
by dividing t by this t-variation. Unfortunately, this variation is 
unknown a priori, so it would be necessary to assume further that the 
Lagrangian scales are directly proportional to the Eulerian scales, about 
which there is previous experimental information. His work has now been 
published) see reference 22. 
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A scale and a microscale can also be defined for R 

/ 


r 



(15) 

(16) 


With nondecaying turbulence, no one-parameter true (time) Lagrangian 
correlation function exists, and the T]- formulation is much more con- 
venient. A further significance of this variable will appear in the 
comparison of Eulerian and Lagrangian treatments of diffusion from a 
line source in flowing turbulence. 


Accelerations in Decaying Turbulence 

A series expansion of v(t + t) for decaying turbulence will show 
something about the initial behavior of the true (time) Lagrangian 
correlation function and will indicate an experimental method for 
examining a hypothesis of Taylor on the interchangeability of instanta- 
neous time and space derivatives when the turbulence level is low 
(reference 23) • 

Write the Lagrangian correlation coefficient 


R v (t, t) 


v(t)y(t - t) 
v'(t)v'(t - t) 


(17) 


Substitute 



into numerator and denominator, and restrict the analysis to small values 
of t : 



(18) 
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Divide numerator and denominator by v, expand the square root in 
the numerator, and keep terms in t^: 


R v (t,r) ~ 1 


dv ' \ ^ + / dv\ ^ T 2 
dt j t \dt/t _ 2 v 2 (t) 


(19) 


For negligible decay rate this reduces to equation (6). 

Equation (19) shows that a Lagrangian microscale defined by 



lim 
t — >0 


1 


Rv(t,T) 


a 2 R v (t,r) 


j n 

1 T =0 


is expressible as 


1 

X ( t ) 




( 20 ) 


(21) 


Introduction of Taylor's ^-postulate transforms equation (19) to 


R T| (n) 


/dv\ 2 _ 

/dv'\ 2 

W n =o 

\ dT ) /n=o 


2v 2 

n=o 


( 22 ) 


after the additional approximation that v'(t)) rs v'(0) when T} is very 
small. 

This gives a new expression for Lagrangian microscale 

.a 


i i 

P r 

^ 2v c 


/dv\ 2 _ /dv'\ ' 
\dri/ “ \dTl / 


( 23 ) 
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Equation (23) is ‘in contradiction to equation (17) of part IV of 
reference 2. In that work Taylor has apparently assumed that v' is a 
nondecaying function of T| . However it certainly is decaying, even in 
terms of this distorted coordinate, and the • dv'/dT)-term must be^- 
included. 

Since and v'(x) can be determined experimentally from the 

mean thermal wake behind a line heat source, equation (23) permits 

/ d \2 

determination of j which is simply related to the mean- square 

"Stokes" acceleration ~ • This quantity is of particular 

interest for the possibility of an instantaneous space-time transforma- 
tion at low turbulence levels. This was first proposed by Taylor (refer- 
ence 23) and has since been used very widely, especially to get approxi- 
mate values of partial derivatives with respect to x (the me an -flow 
direction) by measurements of time partial derivatives. 

The total (or Stokes^ derivative of v(x,y,z,t) in a turbulent 
flow with mean velocity U' along the x-direction is 


dv . dv 
dt dt 


(U + u)|l ♦ 




(2M 


Taylor's hypothesis amounts to the statement that 


( 


dv ^ rr dv 

— ro -U TT 

bt Bx 


or 


bv 

bt 




with 



(24a) 


« 1 


(25) 
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In detail, equation (25) is 


fdv\ 2 6v dv 

V dt) + u i u j a Xi a Xi 


_ dv dv 

2ui — - — 

dt d Xi 


U 2 


1/2 


« 1 


(25a) 


In the absence of information on the algebraic sign of the triple 
correlation term, it is sufficient to require the two conditions^" 


M 


5 "© 


1/2 


« 1 




dv bv 

i J Sxi SxTl 


u,u 


\dx/ 


1/2 


« 1 


(26) 


But 


(S) £ 


can be determined from measurements of Y’Cx)} 


/■\ \ 2 2 

\^-J = 2 where X. is the Eulerian microscale (reference 2)j and 

upper hounds, in terms of measurable functions, can be set on 


u^u^ with the. use of Schwarz inequalities. Thus, an experimental 

J oxi dxj 

check of the requirements in equation (26) is to be made in the section 
entitled "Computation of Results." 

^Lin has discussed the validity of Taylor's hypothesis using a 
slightly different formulation in reference 2b. He points out there that 

n , n 

if A = ^ a a , then (from the Schwarz inequality) ^A 2 j 1/2 ^ n £7) 1/2 

where a^ is a set of numbers. 
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Relation between Eulerian and Lagrangian Microscales 

Taylor (reference 2) inferred an approximate relation between 
and X by neglecting the effect of viscosity on pres sure -gradient 

1/2 

as being approximately 


ratio X^/X for all turbulence. 

The rough nature of this analysis induced Heisenberg (reference 25) to 
conduct a more detailed study of the static-pressure fluctuations and to 
reestimate the X^/X ratio. However, he followed Taylor in ignoring 

the dv'/dT]-term in the relation between X^ and (see equa- 

tion (23)) and in neglecting viscous terms in the relation between X^ 

- w- 

Although these omissions are probably not serious except in the 
low Reynolds number range, it seems interesting, if only for the sake 

of completeness, to use a Heisenberg type of approximation for 

and to repeat his treatment with the omissions rectified. 

From the complete Navier -Stokes equations in the wave-number space, 

n 

Heisenberg deduced an approximate expression for (vp)^ i n terms of 
mean quadruple products of the "harmonics" of the velocity field. His 
principal simplifying assumptions were: 

(a) Different Fourier components of the velocity field are 
uncorrelated 

(b) The turbulent energy spectrum is given by the solution to his 
equilibrium-energy-transfer equation, above a lower cut-off wave num- 
ber k Q 

Following these, but using Chandrasekhar’s (reference 2 6) solution 
to the Heisenberg equation instead of the interpolation formula used by 
Heisenberg, there results 



fluctuations and estimating 

] 1/2 


3p 


f(l)j 


f dp \ 2 


This led to a const 



25.6 

K 



(27) 


where k is a dimensionless empirical constant. 
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The numerical constant in equation (27) would perhaps have been 
given more accurately by the use of a "self -preserving" spectrum 
(calculated by Chandrasekhar from Heisenberg's equation) instead of the 
stationary spectrum with low cut-off wave number. Time was not taken 
to make the requisite additional calculations because: (a) The value 

of (| e) 2 depends principally upon the high-wave -number region of the 

velocity spectrum rather than the low-wave-number region, where the 
difference would be greatest, and (b) the experimental results and 
(especially) the value of k both have a considerable range of 
unc ertainty . 

The mean square of the y-component of the Navier-Stokes equation 
will lead to a relation between and X,: 

dv 
dt 


therefore, 


m - M + ^ 2 *> 2 ■ <*> 

where the correlation between pressure gradient and velocity Laplacian 
function is zero because of isotropy. 

Equations (23) and (27) give the first two terms in terms of the 
microscales, and the mean-square Laplacian function is expressible in 
terms of the fourth derivative of the Von Karman-Howarth f(r) correla- 
tion coefficient at r = 0 (reference 27): 


bv , bv , bv bv 

— + u — + v — + w — 

dt dx dy dz 


1 dp , £ 

■p^ + W 


P oy 


(v 2 v) 2 = ^(v') 2 f iv (0) 


Consequently, equation (28) becomes 



2^.2 

K 


(v') : 

x.3 


32 

3 


v 2(v.)2fiv( 0 ) 


( 29 ) 
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The second term in this equation can be replaced by the turbulence- 
decay equation: 


a(y') 2 

dt 



whence: 


dv' 

dt 



so that: 


(30) 


2 + 25 V 2 ( v, ,> 2 

> 2 +25v “5?T 


^2 + f ^(v')2fiv( 0 ) 

K X.3 3 


therefore. 


2,2 


= 12 , 6 ^ _ 12,5 32 iv( ) 

* R X R X 2 ^ (v*) 2 


(3D 


, r> v'X 
where R x = — . 

Batchelor and Townsend (reference 19) have deduced an expression 
for f iv (0) which is valid in the region of decay where both l/(v') 2 
and \ 2 increase linearly with t (corresponding to large R\) : 


^fiv(o) _ ^0 + i R^S ( 32 ) 

3/2 

y the skewness factor. Their experimental 

results showed S = 0.39 approximately constant for isotropic turbulence. 
Then the estimate for X/A,q becomes 


where S 



m 




(¥♦ 



i_ 

R \ 


( 33 ) 
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The value of k was first estimated by Heisenberg (reference 25), 
from m'easurements of turbulence decay, as O.85. This method may be 
regarded as emphasizing the (relatively low wave number) energy -bearing 
range of the spectrum. Lee (reference 28) worked out an estimate based 
upon skewness factor (k » 0.13), which gives heavy weight to the high- 
wave-number range. Proudman (reference 29) has reestimated k by 
comparison with measured curves of the double and triple velocity correla- 
tions. The value k = 0.45 leads to reasonably good agreement for the 
moderately high -wave-number region, over a wide range of R^. 

It may be remarked that the supposed constancy of k is merely a 
postulate of the Heisenberg dimensional formulation of the spectral 
transfer function. In fact it is by no means obvious that this turbulent 
part of the transfer is quantitatively independent of the amount of 
spectrally local dissipation to heat. In any case, Proudman' s estimate 
of -k = 0.45 has been used here. Therefore, 



( 34 ) 


In the limit of R ^ — > 0, equation (34) does not apply since 
equation (32) does not apply. However, the appropriate limiting rela- 
tion can be obtained directly. In this limiting condition the pressure 
term in equation (28) is negligible compared with the viscous term (the 
former goes with l/R^ and the latter with l/RA 2 )^ and the Eulerian 

velocity correlation coefficient is (reference 27) f(r) = e -r2 / 2 ^ 2 
This gives f iv (0) = 3/A and 



( 35 ) 


Eulerian Analysis of Heat Diffusion from a Line Source 
The two-dimensional turbulent-heat-transfer equation is 


- 80 , ^ 80 k /8 2 0 , 8 2 0\ 
" $ v 57 * + a?) 




( 36 ) 


where 0 is mean temperature, -d is temperature fluctuation about the 
mean, k is thermal conductivity, and Cp is specific heat at constant 
pressure. 
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For the thermal wake behind a line source in isotropic turbulence 
with constant mean velocity V = 0. With restriction to low turbulence 
level, a "boundary -layer" type of approximation can be applied to the 
mean wake, so that 


a 2 © a 2 © 

ax 2 a y 2 


and 


a 

3^ 


(flu) « U ^ 
dx 


so that equation (36) takes the approximate form 


- a© = a3f _ a 

3x pc p ^ y 2 " 3y 


(flv) 


(37) 


It must be emphasized that for this particular initial condition 
on the temperature (effectively a "point source"), the restriction to 
small turbulence level v T /u « 1 does not imply that 3 r /0 is small. 
In fact, for this problem ti'/® is often greater than unity, especially 
at the "edge" of the mean wake, as has been discussed in reference 10 
and will be brought out again later in the present report. 

When the molecular transport can be neglected relative to turbulent 
transport 


u if - - 


(38) 


an equation given in reference 10; a slightly more general treatment 
follows: 

With a constant rate of heat generation (similar to steady state 
in the average), the application of a Von Karman integral-relation 
treatment to equation (38) yields an integral condition: 


2pc p U 



0 dy = Constant = H 


(39) 
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where H is the average time rate at which heat crosses all planes 
perpendicular to U per unit length of heat _source (z-direction) . Of 
course, tfu has been neglected relative to 0U in equation (39). 

Equation (38) has two unknowns, and the first objective is to 
express fir as a function of the (more easily measurable) ©(x,y). 
After integration with respect to y, 


dv = -U 


, y sf dy + F(x) 


But, by symmetry, dv = 0 for y = 0. Hence F(x) = 0 and 


dv = -U dy 

Jo ax 


m 


This relation is sufficient for the computation of dv(x,y) from 
the measured ©(x,y) but the empirical fact of simple geometrical 
similarity in ©(x,y) suggests exploitation of the consequent 
simplification. 

Assume 


©(x,y) = 0 o (x)f(|) (4l) 

where | = y/Y'(x). This transforms equation (39) to 

0 (x)Y'(x) = ~ (42) 

U 


where H* = 


H 


= Constant. It transforms equation (4o) to 


2pc p f(6) d! 

Jo 


dv = [-U Pf(l) dljY'(x)^ + (u P^d^x) 

JO \ Jo 


dY' 

dx 


(^ 3 ) 
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With equation (42), 0 o (x) can be eliminated from equation (43), 
and after integration of the dimensionless integrals this leads to the 
final form for the turbulent -heat -transfer correlation. 


OvU.y) - 


(WO 


where 


t y 

The same sort of analysis can be made on equation ( 37 ) which 
includes the molecular conduction, but the rather large experimental 
scatter in the present measurements seems to make such a refinement 
inappropriate . 

An "exchange" coefficient or "diffusion" coefficient for turbulent 
heat transfer kp is simply expressible in terms of 0(x,y). A con- 
ventional procedure for semiempirical analyses is to write for the 
turbulent transport an expression just like that for the molecular 
transport: 


-pcpOv = |§ 


( 45 ) 


which serves as the definition of krp. 

For the v simple case of equation (MO, it turns out that 


krj. = pc p U 


if(p 

df 

d! 


Y'(x) £ 


(«S) 


This has the particularly interesting property that in a nondecaying 
turbulence at very large values of x, where the mean thermal wake 
spreads parabolically (reference l), kp becomes independent of 

explicit dependence on x. 


A more startling simplification follows for a particular mean 
temperature distribution across the wake: All dependence of km on 


y 
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disappears if 


— a* If 

d| 


that is, if 

_l e 2 

f(|) = e 2 

But this is the Gaussian function, which is found empirically to fit 
all thfe measurements within the experimental scatter. Hence one arrives 
at the empirical result that in both nondecaying and decaying turbulence 
hp is independent of y in the thermal wake behind a line source of 

heat. From equation (4-7), — = -!f(£) and 

kp = pCpUY'(x)^- (48) 

It can be seen that in the nondecaying case at very large values of x, 
kip is constant and independent of both y and x. 


Relation between Some Lagrangian and Eulerian 
Parameters in Transport 

There has apparently been little effort to relate the Eulerian and 
Lagrangian formulations of turbulent diffusion up to the present time. 
Exceedingly simple boundary conditions permit some connection to be 
made in restricted ranges of the present problem. 

For nondecaying or (with much less accuracy) decaying turbulence, 
equation ( 13 ) applies: 


Y 2 = 2 



(R - Rl^nW^i 


fi . ar fd 

dT) dT) 



R R( r >l) dT U 


or, equation (12), 
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For small values of ti. R_ = 1 - -3 — whence 

n V 2 


Y 2 = T) 2 1 


iji 


Si = sJi - 1 

*1 \ 3 


Jill 

X 1 V 


Substituted into equation (44) these give 


(49) 


(50) 


dv = H*f(D^l 


_1_ T) 2 \v» 

12 4 n 2 ju 


(51) 


while equation (46) becomes 


= P c P v ' y 

dl ' 


* 4 ) 


(52) 


“ here ^ ° I Jo v ’ dx “ a 1 = vfa * if + s $)■ In fact > for 

V f 

small values of tu t] a* — x. 

U 

For low-level nondecaying turbulence, T| = ~ x, and equation (51) 
becomes U 


dv = H*f( |)-^ll 


1 v^ _x 2 \ U 
U 2 X r 


(53) 


while equation (52) becomes 


kT = pc p v ' y 


(54) 
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As T| — » 0 both equations (51) and (53) reduce to 


= H*f(S)£ -3L 

X 


while equations (52) and (5^) reduce to 


V - DC V' 

X - PC P V df 

d| 


At the other extreme, when q is very large. 


where M]_ = 


r 

Jo 


= 2T!L n - 2Mi 

tjR^ ( T) ) dr) = Constant and 


In this case, 


^=2L, 

dt) h 


Y L r 


dv = H*f ( | ) — 


U 


f (iW, - Ml)] 


372 


and 


, f (ii 


= pc p v df r 

d£ 


(nLn - M i) 


1/2 


As indicated previously, the ^-variation can be expressed in 
of x or t. 


(55) 


(56) 


(57) 


(58) 


(59) 


(60) 

terms 
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If r| is allowed to become large enough to make yl^ » Mj_ 


v' 


8v = 

(61) 

"T = PC P V ' y (^j 

(62) 


d| 


The above formulas take on particularly simple forms if the 
empirical result of a Gaussian f( |) is utilized: 


fU) 



Then the general expressions for dv and kp (equations (44) and (46)) 
become 


dv(x,y) 


H* dY' 
(Y')2( x ) dx ye 



(63) 


and 


kp(x,y) = pc p v'Y'(x)^- 


(64) 


the latter having been deduced in the previous section. The particular 
forms for small values of t) would follow from substitution of equa- 
tions (49) and (50) into these two. 

However, the most interesting form occurs for very large values 
of q. There Y'(x) is given by equation (57) and 


3v(x,y) = 


H*yLy ^ 




-1 3/2 


exp 


Mn 



(65) 


kp(x,y) = PCpV'Lri - Constant 


( 66 ) 



34 


NACA TN 2710 


For still larger values of 


tj, such that M]_ « t^L^, 


and 


3v(x,y) 


y xl 

H* * U 

exP 



(67) 


fc|,(x,y) = pCpV'L-ri = Constant 


The constancy of kip for large values of r| ( large values of t 
or x) is to be expected} a treatment of molecular diffusion by this 
method must certainly yield a constant coefficient for times much larger 
than the mean free time (of flight) of the molecules - that is, for all 
"macroscopic times." Put another way, the simple parabolic_behavior 
of Y' for large values of T} is a sure indication that 0 obeys the 
simple classical diffusion equation with constant coefficient, when 
viewed extremely "coarsely." 


Perhaps the chief 'interest of equation (66) is its identification 
of Lq as a significant Lagrangian length for diffusion at a large 

distance from the source. It enters the expression for turbulent- 
diffusion coefficient in much the same way as mean free path enters 
the expressions for the molecular-diffusion coefficients. Furthermore, 
its role appears to be much like that attributed to Prandtl's "mixing 
length," which was brought into the turbulent -transport problem in a 
more or less intuitive fashion. 


Of course, the possible crude nature of Taylor's original 
^-postulate may render the significance of Lrj more qualitative than 
quantitative in the case of decaying turbulence. 


COMPUTATION OF RESULTS 


Although Taylor's assumption of the unique dependence of t^t-T 
upon T| is not likely to be accurate for collapsing together cases 
with widely differing turbulence decay rates, it does provide a rela- 
tively simple relation between Y'(x) and Rt)(r). Therefore all of 
the mean- thermal-wake data were reduced on the T]-basis. 

In principle the complete R^(t]) curve can be obtained from Y^(x) 
by double differentiation (equation (12)): 
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or 



( 68 ) 



However, simple double differentiation of the squares of a curve 
as uncertain as Y'(x) seems almost hopelessly indeterminate - altho ugh 
Taylor (reference 2) and Collis (reference 7) have apparently followed 
this procedure. A somewhat more circumspect technique has been tried 
here: The values of and were determined first, thro ug h certain 

limit relations (to be described) . Then the R^ curve was determined 
by double differentiation, subject to the restrictions of agreement with 
the previously determined scales. Thanks to rather poor determinacy of 
values of A-t) and this method is not so much of an improvement as 

it might first appear. 


Lagrangian Microscale Xt) 

If equation (13) is restricted to very small values of q the 
parabolic approximation for R^ can be introduced: 


Y 2 = 2 





therefore 


(69) 


whenc e 



£ jL 

<^ 2 



(70) 


The computational procedure was to plot 
to estimate the slope of the faired curve at 


Y 2 /t) 2 against t] 2 and 
T] = 0 where the curve 
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must pass through unity. The abscissa intercept of the O-tangent is 
6X T j 2 . The actual points, faired curves, and tangents for all cases are 
presented in figures 13 to 22. Clearly the precision is poor. 


Lagrangian Scale 

Consider equation (12) in the limit as T| ^oo. It imm ediately 

gives 

V-5 (if) (7D 

and the graphical procedure based on this is also presented in figures 13 
to 22. Some of the asymptotic slopes drawn are not the best representa- 
tion of the experimental points. This is due to the auxiliary (assumed) 
restriction that Rtj cannot increase with increasing values of tj as 
long as R^ has not previously dropped below zero. The graphical 
precision attainable is perhaps a little better here than that for X^ 2 , 

but the square root necessary to get X^ means that X^ is determined 
about as well as is 1^. 


Lagrangian Correlation Function R^(t]) 

With X-q and determined, the initial (small Ax) and final 

(large Ax) behavior of the curve Y’(Ax) is prescribed. These parts 
of the curve were drawn on a graph with the experimental points. Then 
the fairing in of a reasonable central portion to this mean Y'(Ax) 
curve was a relatively simple matter. The Rrj curve was then obtained 
by double differentiation. 

The curves drawn for Y'(Ax) in figures 23 to 32 were determined 
in the fashion described above, as were the curv.es for R^ in the same 
figures. 


Eulerian Microscale X 

In view of the approximate nature of the determination of X^, no 
new direct measurements were made of the Eulerian microscale X. Instead, 
X was computed with the energy equation for isotropic turbulence, from 
the measurements of turbulence decay:. 
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2 


-10 V 


dv 2 

dt 


(72) 


or, with the space-time transformation. 


X = -lOv 


v 2 


(73) 


Eulerian Scale L 

Earlier investigations have shown that the Eulerian scale in a 
grid-produced turbulence is closely a linear function of the mesh size 
of the grid producing the turbulence (for. a given value of x and grid 
geometry) and is not significantly dependent upon the mean velocity 
(or grid Reynolds number, provided it is sufficient to cause turbulence). 
Therefore the values of L have been deduced from earlier measurements 
at the California Institute of Technology (reference 18) on grids of 
essentially the same geometry. 

Table 1 summarizes the results for Lagrangian and Eulerian scales 
and microscales. The results have been grouped to show the effect of 
systematic variation of one parameter at a time. Some of the results 
are presented in figures 33> 3^.> and 35- 


Instantaneous Space-Time Transformation 

The permissibility of an instantaneous space-time transformation 
in flowing turbulence. 
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can be estimated in accordance with equations (26). For equation (7^) 
to be valid, the sufficient requirements are those given in equations (26) 
that: 




With the aid of equation (23) 

/ 6 v\^ v‘- 

the Taylor relation = 2 — , 


, the turbulence decay equation, and 
the first of these conditions can 


be written in the form 


M 

12-5 + M 2 

W 

% 2 W 


1/2 


« 1 


( 75 ) 


From 

magnitude 


the Schwarz inequality, essentially the necessity that the 
of any correlation coefficient be less than or equal to unity. 


dv dv < 

u i u j sq; zr; “ “i 


•u. 




( 76 ) 


where the prime in this expression denotes root -mean- square value. For 
isotropic turbulence, equation (76) can be written 


u i u j 


dv dv < ( n 

^7^7= l 9 


H 2) 


(v ) 


( 77 ) 
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Thus, the second condition in equation (26). will he satisfied if 

2.7 y«i m 

Both T and v'/U for the flows studied are presented in table 1 . 
It is clear that for these flows instantaneous x and t partial 
derivatives may be taken proportional with reasonable confidence. 


dv-Correlation 

The Eulerian measure of transverse turbulent heat transport is 
computed from the mean temperature distribution. The dimensionless 
form, i 8 v/© 0 U, is given for two typical cross sections in figures 36 - 
and 37 - 

The measurements of v'/u and of d '/© 0 permit calculation of 

v 

the correlation coefficient and this is also given in 

■Q V 

figures 36 and 37 * 

For the convenient and reasonably accurate assumption of Gaussian 
mean temperature distribution, the corresponding turbulent-heat-transfer 
coefficient bp follows from equation ( 6 k). It was found to be inde- 
pendent of y, and typical curves of bp/k are given in figure 38 . 

The data for bp at three different speeds behind the 1-inch grid are 
roughly collapsed together through division of kip by pcpV'L^, as 
suggested by equation ( 66 a), an asymptotic result for nondecaying 
turbulence (fig. 39 ) • 


DISCUSSION 
Lagrangian Variables 


Even a cursory examination of the technique used in this investi- 
gation for the determination of Lagrangian correlation shows that, as 
physical measurements go, this method is a "bad" one, largely because 
of the inherent double differentiation between measured variable and 
desired information. 

Figures 13 to 22 suggest an uncertainty in values of and 

as large as ±20 percent, in spite of moderately good precision in the 
measurement of individual temperature distributions such as the lower 
curve in figure 3 - 
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As mentioned earlier, the values of v'/U computed from initial 
wake spread are consistently higher than those measured with the hot- 
wire anemometer (reference l8) . The same relative result was encountered 
during a brief investigation following that reported in reference 18 . 

Up to the present time there has been no satisfactory explanation of the 
discrepancy. A tentative hypothesis which would at least account for 
its direction may be based upon a human weakness in the visual averaging 
of the reading of a fluctuation pointer^ there seems to be a tendency to 
choose an "average" more or less halfway between the extremes of the 
needle travel. Thus a pointer motion with very skew probability density 
(greater than 0) would tend to be "averaged" at too high a value. The 
thermocouple voltage in one of these thermal wake traverses has just this 
character (fig. 11) . Hence a visual averaging might yield too high a wake 
width. If this effect is nonnegligible, it is advisable to employ some 
electrical means of averaging for skew signals, for example, the flux 
meter and bucking circuit described in reference 10. 

In view of the considerable uncertainty in as well as that 

in X and v* the poor degree of agreement between experiment and theory 
shown in figure' 33 is understandable. Since the two undetermined con- 
stants in the theoretical result have been evaluated from sets of experi- 
ments completely independent of the present ones, this agreement can be 
viewed as an affirmative result. 

Since some sort of Lagrangian scale should be a significant length 
in turbulent heat and mass transport, as demonstrated in the analytical 
section of this report, an effort has been made to find some systematic 
variation in the values of L^. Figure 3^- might be construed to indicate 
a monotonic decrease of Ltj/L with increasing values of R^. It is 
interesting to note that a decrease was also observed for the ratio of 
mixing length to tube radius by Nikuradse (reference 30) in fully devel- 
oped turbulent tube flows. In order to determine whether these two rates 
of decrease with increasing Reynolds numbers are of the same order of 
magnitude, an estimate has been made of the magnitudes of Rl corre- 
sponding to Nikuradse' s results given in figures 28 and 29 of refer- 
ence 30. Both scale-to-diameter ratio and average turbulent levels for 
various Reynolds numbers were estimated with the help of Laufer's data 
on turbulent channel flow (reference 31) at various Reynolds numbers. 

The absolute level (i.e., the ordinate scale) of the resulting Z mav /L 
against R^ curve was adjusted to give the most reasonable-looking fit 
with the L^/L data. This is the dashed line in figure 3^ and it shows 
at least a qualitative resemblance. It is likely, however, that L^/L 
is not a unique function of R^. 

Of course not all of the scatter in figures 3^- and 35 (involving Lr^) 
can be attributed to simple lack of experimental precision. Some is 
evidence of the fact that the Taylor postulate of Lagrangian correlation 
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function being uniquely a function of t} is certainly not very closely 
true. Furthermore, table 1 does show rather systematic variations of 
in some of the three-point groups. Most noticeably, there is a 

regular decrease in with increasing x 0 /M (or perhaps with 

decreasing v'/u) for each of the three grids. 


Temperature-Fluctuation Field 

Fairly cl^se behind a line heat source in turbulent flow, the 
random pulse nature of the temperature fluctuations at a fixed point 
has been established by the oscillogram in reference 10. This is 
confirmed by the first two oscillograms in figure 11, with their highly 
skew probability densities at Ax/M = 10. One of the objectives of the 
present investigation was to find out whether this distinctly pulsed 
character persisted far downstream or whether molecular heat conduction 
becomes increasingly effective in smearing out the pulses, until they 
are no longer distinguishable as such. The third oscillogram and proba- 
bility density in figure 11 (Ax/M = 70) does show a decided trend away 
from the pulse-type signal. The molecular broadening of the laminar wake 
(corresponding to the pulses) decreases the relative spacing of the pulses 
•in d(t) at any point in the "turbulent wake" region. This is a reduc- 
tion in relative length of the flat (d = 0) base lines between pulses, 
giving greater statistical symmetry in d(t) about its mean, that is, 
reducing the skewness of P(d). 

A simple analysis will show the existence of an asymptotic behavior 
of molecular-conduction effects in a nondecaying turbulence.. For a 
nondecaying flow turbulence and very large values of t(~ x/U) the mean- 
square wake spread due to turbulent motion is 


Y 2 ^2L t1 ^-x (79) 

u 

On the other hand, an approach to molecular diffusion through 
Taylor's concept of "continuous movements" gives, for any macroscopic 
distance downstream, the mean-square thermal wake width, 



U 


(80) 


where A is the mean free path and c -is the root -mean- square molecular 
velocity. 
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From equations (79) and (80) 





(81) 


For a typical case ; take L = 1 centimeter, v f = 10 centimeters per 
second, A = 6 X 10“^ centimeter, and c = 5 X lcA centimeters per 

l — ~ 0.17. For people accustomed to thfnking of molecular 


second. Then 


transport as negligibly small in turbulent flow (e.g., in shear flow), 
this ratio will appear quite large. The values of kp/k plotted in 
figure 38 also show that at these low turbulence levels the molecular 
thermal conductivity is not necessarily negligible compared with the 
turbulent transport . 


The temperature-fluctuation-level distributions across the 

wake (figs. 10(a) and 10(b)) show the same character as that measured at 
much higher turbulence level in a jet (reference 10), with somewhat 
lower minimum values, which are attributable to the lower turbulence 
level. A rough evaluation of the behavior of the statistical variables 
in this turbulent thermal wake is obtainable by recalling that it consists 
of a randomly "waving" laminar thermal wake. If $(t) is crudely 
represented by a randomly spaced sequence of identical rectangular pulses 
with height $ 0 , width j, and average spacing s, it is easily seen that 



This permits -a 1 /© to vary between 0 and 00 as s/j travels the 
permissible range from 1 to 00 . Since points nearer the edge of the 
turbulent thermal wake have higher values of s/j, the behavior of 
.equation (83) is consistent with the experimental distribution. If 
the analysis were repeated with triangular pulses, for example, the 
quantitative estimate would doubtless be more realistic. The higher 
values of ('3 T /50min encountered at higher values of v T /u are an 
indication that for a given width of laminar thermal wake, the higher 
value of v ! /u leads to a higher minimum value of s/j. 
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The closely Gaussian shape of 0/0 o against y has already been 
pointed out . If 


#- = ex ? 
u o 


. i 

2 (yf 


be introduced, there results 


sL =< 
8 ' 

do 

— exp 
©0 * 

1 y 2 . 

_ 2 (yf 

- ll 


L 




1/2 


and 


m 


d' 

nr exp 


i jd_ 

2 (vf 


d 0 

— exp 
0 


( 


i jd. 

2 W 


- 1> 


1/2 


(85) 


Both of these expressions have behavior consistent with the experiments. 


The form of dimensionless transverse turbulent -heat -transfer rate 
dv/0 Q U can be deduced for small values of Ax (such that ~ 1) with 
this pulse representation of d(t). In this picture, fv is the correla- 
tion between a continuous random variable v and a random pulse signal 
d which "fires" every time the continuous variable passes through a 
specific value 


Therefore 


■t-X = jr~ U 
1 Ax 


dv 


dp J 
s 



_U_ 

Ax 


oy 


( 86 ) 
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where s and 0 are functions of y. 


With Gaussian 0(y), v 


$ v _ y_ 

ioU _AX 


exp 



(87) 


The direct comparison between this crude picture and the experimental 
results will be confined to the correlation coefficient R^ v = ^v/d'v. 
This is of particular interest in view of the surprisingly high experi- 
mental values. With equation ( 85 ) and the fact that v'/u = Y'/Ax, 
there results 


r 



This contains the undetermined constant d o /0 o , which can be obtained 
from any one of several experimental results. Figure 36 includes one 
plot of equation (88) with d o /0 o determined from equation (84) and the 
experimental value of (d ' /<5) m ^ n , an< ^ one Pl°"k with $ 0 /q determined by 
matching equation (88) with the experimental result at y/Y' = 1.4. 

It is also surprising to find that the experimental (^dv) ffl ax 

large values of Ax is even larger than that at small values of Ax. 

This may be due to a considerable experimental error; the resistance- 
thermometer voltage signal is much lower here. Unfortunately no relation 
corresponding to equation '(88) has been deduced for large values of Ax/ 
where R^ is essentially zero. 


The criteria for Taylor's hypothesis of the interchangeability of 
instantaneous space U ^ and time. ^ derivatives (assumed by him to 

depend only upon turbulence level) have been expressed in equations (75) 
and ( 78 ) as functions of turbulence level R\ and X/X^. If X/X-^ is ■ 
replaced by its theoretical expression (equation (34)) in terms of R 
equation (75) becomes 





(89) 


For estimates of most flows the first term in the parentheses can be 
neglected; values of R^, less than 5 or 10 are rare. 
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Since equation (3^-) has now been roughly verified by experiment, 
equations ( 89 ) and ( 78 ) may serve as approximate criteria for the 
validity of Taylor's hypothesis. 

In the limit of R^ — >0 when equation (35) replaces equation (3*0* 
there follows a simpler criterion to replace equation ( 75 ): 

Ifc-) « 1 (90) 

R x\u/ 


SUMMARY OF RESULTS 


The following results were obtained from the investigation of the 
diffusion of heat from a line source in isotropic turbulence. 

1. The thermal wake behind a heated wire set perpendicular to a 
flowing isotropic turbulence (at sufficiently low wire Reynolds number) 
consists of a randomly "waving" thin, laminar, thermal wake whose varia- 
tions in lateral position give what may be called the turbulent thermal 
wake. At a ^fixed point not too far behind the wire t 3 ( t ) is a random 
pulse function, and' the nature of the turbulent heat transfer can be 
deduced on this basis. Farther downstream the distinct pulse nature 
tends to disappear. 

2. The mean transverse temperature distribution -JKy) appears to 

®o 

be Gaussian within the experimental precision for all distances behind 
the wire. 

3. An Eulerian analysis of this turbulent -heat -transfer problem 
permits computation of the turbulent -heat -transfer coefficient kp < 
which is essentially constant with respect to the distance in the 
direction of the measured diffusion y for these boundary conditions. 

It is found that at low turbulence levels (approximately equal to 1 
to 2 percent) the molecular heat transport is not vanishingly small 
compared with the turbulent heat transport. 

h. Although Taylor T s postulate that Lagrangian correlations in 
decaying turbulence can be made similar by introduction of an independent 

ft 

variable r\ = v*(t) dt (where t is time and v* is the root- 

mean-square instantaneous velocity fluctuation in the y-direction) 
seems to be an oversimplification, it has been applied here for 
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convenience in the reduction of data. A simple comparison of Eulerian 
and Lagrangian analyses for diffusion in nondecaying turbulence shows 
that for large values of the distance from the heat source Ax the 
Lagrangian scale enters the expression for kip, the turbulent- 

heat-transfer coefficient, much like the empirical mixing length in 
the old turbulent transport theories. Therefore some properly modified 
generalization of Taylor's T] -postulate should prove useful. 


5- A correction and generalization of Heisenberg's theoretical 
expression for the ratio of Eulerian to Lagrangian microscale X/Xf) 
as a function only of the turbulence Reynolds number based on microscale 
R\ has been made and seems to agree roughly with experiment. It must 
be noted that since X^ depends only upon a transformation dr) = v' dt, 

and not upon the integral postulate stated above, its validity is not 
impaired by any failure of the integral postulate. 


6. Taylor's hypothesis for the interchangeability of space and 
time derivatives at low turbulence levels has been expressed in terms 
of criteria which depend upon turbulence level, Reynolds number, and 
X/Xt). Applied to the flows studied here it shows that in these cases 
such a transformation is permissible. By substitution of the theoret- 


ical expression for ^-(r^), a slightly simpler and rougher criterion 


is derived, depending only upon turbulence level and R-\. 


The Johns Hopkins University 

Baltimore,, Md., June 5, 1951 
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(a) No grid in tunnel. 



Figure 2.- 


(b) Grid-produced turbulence. 



Shadowgraph time exposure of wire wake. 




38 ft/sec 
25.6 ft/sec 



Figure 4.- Spread of heat from a line source. — = 43.4, M = 1 inch 
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Figure 9.- Temperature fluctuations at x Q = >3 inches, M = 1 inch, and 
U = 25.6 feet per second along wake axis. 
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(a) Spread of heat from a line source. 
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(h) Correlation function R^. 


Figure 29.- Spread of heat from a line source and. correlation function R^ 

1 — x o 

for M = — inch, U = 25.6 feet per second, and — = 86.8. 

2 M 
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(b) Correlation function R^. 

Figure 30.- Spread of heat from a line source and correlation function R^ 

for M = i inch. U = 25.6 feet per second, and — = 86.8. 

^ M 









Figure 33-- Ratio of Lagrangian and Eulerian microscales 
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Figure 35.- Ratio of Lagrangian microscale and scale. 
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(,aj — — - against — , 

e n u y* 


^ — ^Free constant determined 

v Free constant N ' . 

y ^ 

determined from R. at A=|.4 


(b) R^ v against y/Y*. Dashed curves are result of analysis based on 

rectangular temperature pulses. 

Figure 36.- Heat-transfer correlation across_thermal wake, computed from 
measured mean temperature distribution. U = 25.6 feet per second, 

• M = 1 inch, x Q = 43.4 inches, and Ax = 10. 5 inches. 












Figure 39.- Dimensionless turbulent-heat-transfer coefficient. 




